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of the scalar field and gravitational Lagrangians. 
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1. Introduction 

The number of spatial dimensions of our Universe is a key foundational question that 
has received enormous theoretical and experimental attention in the last decades [1]. 
Initially inspired by developments in string theories, the idea of braneworlds as four¬ 
dimensional hypersurfaces embedded in a higher-dimensional (typically 5—dimensional) 
bulk space has been explored with exquisite detail in the literature Among the 

resulting constructions we find particularly interesting the case of thick branes. This 
setup consists on a bD bulk space with a scalar held localized in the hfth dimension 
and whose energy density is typically concentrated around the point ?/ = 0 of the extra 
spatial dimension, where the 4D brane resides. The rest of the matter helds are conhned 
on this (thick) hypersurface, with the scalar held somehow acting as the force responsible 
for their conhnement. This provides an interesting framework, where both the fermion 
mass hierarchy and proton stability can be investigated due to the thickness of the 
brane; see, e.g.. Ref. |1]. 

The notion of a higher dimensional space-time is so fundamental that transcends 
its original motivations and justihes its exploration from diherent perspectives. We 
thus hnd it useful to seek for inspiration in a lower dimensional analogy borrowed 
from condensed matter systems and interpret braneworld models as a kind of sandwich 
conhguration in a crystal, with a (possibly thick) 2—dimensional distribution of 
impurities within a 3—dimensional crystalline structure, being the crystal analogous 
to the bulk and the (thick) layer of impurities the brane. Interestingly, the structure 
of crystals in the continuum limit admits an effective geometric description mum- 
This point is relevant in our discussion of braneworlds because the presence of defects 
such as interstitials (either impurities or not) or vacancies in crystals requires the use of 
metric-affine geometry for a consistent description of their effective geometry. Indeed, 
a direct correspondence between the metric-affine geometry of crystals with point-like 
defects (interstitials and/or vacancies) and theories of gravity with non-metricity has 
been recently established jH]. It turns out that the density of point-like defects in a 
crystal is analogous to the density of energy-momentum in the gravitational context, 
both being responsible for the existence of non-metricity, a property associated with the 
afhne connection. The conhnement of energy-momentum at or around a region, dehning 
in this way the brane, is thus analogous to the concentration of point-like defects on 
a sandwich layer in a crystal. This analogy, therefore, provides a novel motivation for 
the study of braneworlds in geometric scenarios with independent metric and afhne 
degrees of freedom. Given that the question of whether the space-time geometry is 
Riemannian or otherwise must be determined by experiments, rather than by tradition 
or convention, we believe that the impact that metric-afhne geometry could have on 
the phenomenology of higher-dimensional models of the Universe should be explored in 
some detail. This is the main motivation for this work. 

The dynamical laws governing the higher-dimensional scenario we are about to 
consider are, in principle, unknown. Nevertheless, it is reasonable to expect them to be 
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close to those of General Relativity (GR). For this reason, and to allow for some degree 
of generality, we may assume that the gravity Lagrangian is some function of the space- 
time metric and the Ricci tensor, F^gap, Rag), being GR the case F{gag, Rag) = g^^Rag- 
We will see that, unlike in the (more standard) metric formulation of such theories, the 
held equations for this type of theories are always second-order and ghost-free, which 
is an important feature for quantum investigations. This fact allows us to explore in 
detail the equations governing the evolution of tensorial perturbations on the brane. We 
hnd that with an appropriate choice of variables and some standard redehnitions, these 
equations can be put in a form which is formally identical to that found in the case of GR. 
Their stability is proved by showing the positive dehniteness of the effective Hamiltonian 
associated to the higher-dimensional modes. This result conhrms the robustness of the 
braneworld scenario against tensorial perturbations for arbitrary gravity Lagrangian of 
the form F{gag, Rag) and for arbitrary scalar held dynamics. 


2. Dynamical content 


We consider a rather general family of theories of gravity based on the metric and an 
independent connection (metric-affine or Palatini approach [9]), with action 


S = 


1 


d^Xy/^F{gag, Rag) + , 


( 1 ) 


where S^j) represents the scalar held sector, is a constant with suitable dimensions, D is 
the number of space-time dimensions and g is the determinant of the space-time metric 
g^y. The connection appears through the Ricci tensor, which is dehned as Rg^ = Rt^gau, 
where R°‘g^^iu = dif^vg — Given that the Riemann tensor 

is antisymmetric only in its two last indices, it is not a priori guaranteed that for an 
arbitrary connection T"^ the Ricci tensor has any specihc symmetry. For this reason 
one must bear in mind that there might be an antisymmetric part in Rp^. Variation of 
the above action leads to 


6S = 


1 


d^Xy/^ 






dF 

dRag 



+ , ( 2 ) 


where the variation of the Ricci tensor is given by 


SRg, = (3) 

and = (rj„ - rtg/2 is the torsion tensor. In our metric-affine scenario, the held 
equations follow from equating to zero the independent variations with respect to the 
metric and the connection. Given that 6Rag only depends on the connection variation, 
the metric held equation can be readily extracted from ([2]), which yields 


dF 





( 4 ) 
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where the scalar field stress-energy tensor has been included. The 

connection equation is obtained by inserting ([3]) into ([2]) and integrating by parts the 
first two terms. After elementary manipulations, we find 


Va = 2 ^ 




+ - 


(D-l) 




'ap' 


(5) 


where = dF/dRpi,. At this point one must decide whether to keep the torsion 
and continue with the analysis of the field equations in full generality or introduce 
some simplification to consider specific cases of interest. In this sense, the analogy with 
condensed matter systems tells us that torsion is the representation in the continuum 
limit of line-like defects such as dislocations and disclinations 0 m- If torsion is 
kept in the field equations, the resulting branes could be seen as entities engendering 
two or more extra dimensions able to accommodate line-like defects. Since we are 
focusing on a single extra dimension [21 [10], we find it appropriate to deal only with 
the simplest kind of defects, point-like, which are associated with the non-metricity 
properties of the connection. We leave the consideration of torsional effects for the 
future, where a natural direction should concern braneworld scenarios endowed with two 
extra dimensions mm In the case of a single extra dimension, the above equation 
becomes 


Vj = 0. 


( 6 ) 


A particularly interesting solution of this equation concerns being a symmetric 

tensor, which implies that Rpi, = R^p. In that case, as shown below, the connection can 
be immediately solved as the Christoffel symbols of an auxiliary metric which is related 
to the space-time metric g^^i, through a transformation that depends on the stress-energy 
tensor of the scalar field. To see this explicitly, consider that the Lagrangian F{ga/ 3 , Rag) 
is written in terms of the matrix = g^^R\y and traces of its powers. This type 
of Lagrangians describe the case of f{R) theories [13], where R is the trace of M, 
f{R,R^yR^'^) theories [E], where R^yR^'^ is the trace of M^, Born-Infeld like models 
[laiisi, and other theories considered in the literature. Using the matrix M, we find 
that 


dF 

dgPU 

dF 

dRpy 


1 
2 L 


{Fm) RyX + {Fm)u^RpX 


9^\Fm)x^ , 


(7) 

( 8 ) 


where {Fm)x = ■ In some cases, the symmetry of Rpy follows directly from the 

torsionless condition (see ini). though in general it should be seen as an additional 
simplification, which we assume from now on. Using (I7|) and the definition of M to 
replace Rpy by g^xM^y, the metric field equations can be written as 


{FM)p^gunM\ + {FM)y^gp.M\ 


:9pi 


,F = . 




(9) 


These equations allow to obtain M as an algebraic function of Tpt'^ and g^y. As a 
result, the term ^/—gW^’^' in (jHI) does not depend explicitly on the connection, which 
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allows to find a solution for using elementary algebraic manipulations. In fact, the 
identification 

V^g^\FM)x = , (10) 

leads to 

Aa. (H) 

where \Fm\ represents the determinant of the (invertible) matrix {Fm)\^- This turns 
the connection equation (j6]) into the well-known form 

Va [V^q^’^] = 0 , (12) 

which implies that the solution for T^j, is given by the Christoffel symbols of the auxiliary 
metric This result also means that = R^y{q). Note that via Eq. flTT]) the relation 
between g^y and q^y is purely algebraic and depends only on the scalar held stress-energy 
tensor. 

When T^t'^ is diagonal, the metric gfj_y can be written in diagonal form and (1^ 
suggests that M and Fm also inherit that structure. One can then write (|9]) in the more 
compact form 

= + . (13) 

By raising one index on this equation with and using flTTll . we dually get 

R."{q) = + tW“) , (14) 

\Fm\^-^ 

where Cq represents the gravity Lagrangian of the theory under consideration. This 
Einstein-like representation of the held equations is valid for GR, f{R), and many other 
theories [HI [H [ig El]. It puts forward that qap satishes second-order equations and, 
given the algebraic relations fllip . it follows that the dynamics of gag is also second- 
order. When Ty'^^'^ = 0, fET|) recovers Einstein’s equations in vacuum (with possibly a 
cosmological constant term), which are clearly ghost-free. This last point can be verihed 
by setting to zero, raising one index in Eq.Q with the metric g°‘^, and computing 
traces of powers of the resulting object. Those traces represent algebraic equations 
that relate the scalars R = R^yR^^ = and so on. With a sufficient 

number of equations, which depends on the space-time dimension, one concludes that 
all such scalars in vacuum must be constants. As a result, the matrix M°‘g must be 
proportional to the identity matrix, making in this way = ag^y, with a being an 

(irrelevant) constant. One then readily verihes that Eq. fETll coincides with Einstein’s 
equations in vacuum. 
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We now assume the space-time to be (d-l- l)-dimensional, with a single extra dimension, 
and write the line element for the metric in the form [ 2 ] 

ds^ = af{y)r]abdx°'dx^ -|- dy^, (15) 

being r]ab the metric of a d—dimensional space of constant curvature K. In a thick brane 
scenario, one assumes that a scalar held lives in the extra dimension with most of its 
energy density conhned around the hypersurface |/ = 0, which dehnes the (thick) brane 
[iniiiHiiin]. Here the Lagrangian density in the scalar held action, = J d'^^^XyJ—gC, 
is supposed to be a function C = £(0,X), with X = g'^^dafdpf). This form allows 
that the scalar held engenders generalized kinematics, including the case of k-helds I2D1. 
Given that 0 = 0(|/), one hnds that the stress-energy tensor of the scalar held can 
generically be written as 

— i '^+iy)^dxd G 

0 T_{y) J’ ^ ^ 

where Idxd is the d x d identity matrix, T+ = —£(0, X)/2, and T_ = Cx4>l + d+, with 

the dehnition Cx = dC/dX. As explained above, the diagonal character of induces 

a diagonal which allows us to dehne 

Sjy ^ f f j , (17) 

with and ^l± being functions of (j) and X. We can thus dehne an auxiliary 

line element for of the form 

ds^ = d^{y)r]abdx°'dx^ + dy"^, (18) 

where d‘^{y) = a‘^{y)Q^ and dy"^ = Q_dy‘^ follow immediately from Eq. fll7p . Using 
this to construct R,y^{q), we hnd that RJ' = [{d — 1)K — (Hy + d ■ H'^)]6j’ and 
Ry^ = —d ■ {Hy + H^), where H = a^/a, and K is the (constant) curvature of the 
d-dimensional brane. Inserting this in flTT)) . the background equations can thus be 
written as 


d{d-l)[K-H^] 

{d-l)[K + Hy] 


|G|V2 


{d - 1 )Cg + d ■ T+ - T_ 


K 


|G|V2 


(^4 


T_) . 


(19) 

( 20 ) 


A hrst-order formulation is possible for the case X = 0 if we dehne a superpotential 
W{(j)) such that H = —W{(j))/{d — 1) [21]. The case X 7 ^ 0 is more complicated and 
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we put it aside, recalling that it could perhaps follow the lines of [22]. For 77 = 0 the 
above equations become 


Rx4>y = 


0'^/2 a 




|ff|V 2 


Cxfl+{d-l) 


£( 0 ,X) 


-c 


G 


.( 21 ) 


It is easy to see that for the Einstein-Hilbert Lagrangian, Cg = R/2k‘^, the above 
equations recover the well-known results of GR (in H = d -|- 1 = 5). In this case, 
we get = 1 and Cq = — 5£/2)/3, which lead to Cx4>y = and 

|hF^(0) = K^{2Cxf'l—R)- For the particular case of a canonical scalar held, C = X+2V, 
from this last relation we get 2K^V{(f)) = — |hF^, which is in agreement with the 

results presented in [181CH] . In the more general case in which the scalar held Lagrangian 
can be split into a kinetic part plus a potential term, the hrst of the above equations 
can, in principle, be used to isolate V (0) as a function of (fy and W^. Inserting the result 
in the second, we hnd an (implicit) hrst-order equation for cfyy as a function of IF(0) 
and W^. The fact that the scalar held is generically governed by a hrst-order equation 
suggests that supersymmetric extensions could exist for this type of theories. 


4. Tensorial perturbations 

Using Gaussian normal coordinates, the ehect of tensorial perturbations on the line 
element flT^ can be parameterized as 

ds^ = of{y) {r]ab + hat) dx^^dx'" -7 dy'^ , (22) 

where 5gab = a^{y)hab and 5gay = 0 = 5gyy account for the tensorial perturbations. For 
the auxiliary metric, we must thus have 

ds^ = df{y) {rjab + hab) dx°‘dx^ -|- dy^ . (23) 

Since the tensorial perturbations are restricted to the brane modes, any contraction of 
the form will be zero, as can be easily understood from the structure of (|22|) 

and the fact that only dycj) is nonzero. As a result, perturbation of flT4l) leads simply to 

dR/(g) = 0 O 6R,,{q) = R/tg, , (24) 

where tab = of hab is the only nonzero component of tg^. Using standard covariant 
perturbation methods, one hnds that 


dRfj,u(^q^ 


- \ (V/.VAt^ + V^VAtJ) 

— (P +“ _l_ F? bd 


(25) 


where t^ = t^ — As usual, we will adopt a traceless and transverse gauge, where 
f" = 0 and Va^^ = 0. The transverse gauge can also be written as = 0, with 
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Vq referred to the metric qab, and as Vjh* = 0, with Vj referred to the metric rjij and 
h) = hikv’^R 

From the diagonal strnctnre of (lT6l) and Eq. (IHl), it follows that can 

be written in a diagonal form analogons to ffT6|) and mill- Together with the fact 
that t^ only has non-zero contribntions in the brane sector, one can easily verify 
that I {RkLotf, + Ruatf) = Rfftpy. As a result, the equations governing the tensorial 
perturbations boil down to (we note that similar manipulations can be used to study 
tensorial peturbations in cosmological backgrounds |23]) 


= 0 . 

(26) 

Remarkably, this result is valid for any of the form (lUit. 

lengthy) calculations, we hnd 

With elementary (though 

^^^^'>R^anb = ^^^R^anb + H^S^qan " C^afc) 

(27) 

=dyytl 

(28) 

q^^VmVntl = + dHdytl - 2HX , 

(29) 

where H = hy/a. With these results, fl2^ becomes 


a„ti + dHditi +<“' Dti + = 0 . 

(30) 


It is customary in the literature to model the brane as a maximally symmetric space of 
constant curvature K = ±1,0, in analogy with the properties of our universe at large 
scales. Following this idealization, we get ^'^^R^anc = —K{5f'rian — S^rjac), which turns 
the last term of the above equation into —2Kt^^/of. It is also common in the literature 


to introduce a new coordinate such that the line element ( 12 ^ becomes 

ds^ = of \{'qab + hab) dx°‘dx’^ ± dz^'\ . (31) 

By doing this, fl30l) takes the hnal form 

+{d- l)Hd,ha^ - 2Kha^ = 0 , (32) 

where is computed using the rjij metric and Ri = a^/a. This last equation is formally 
identical to that given, for instance, in [ 2 l] for the case of GR in 4 -|- 1 dimensions. 
Assuming that ha = X{z)ecf‘{t,x), (13^ splits in two equations of the form 

+ mX^ + p^X = 0 (33) 

= 0 , (34) 


where p^ is a constant. At this point, it is important to note that the ea^{t,x) part 
of the tensorial modes satishes an equation which is identical to that found in GR. 
This equation only depends on the brane coordinates {t, x) and the brane metric 
rjab but is insensitive to the gravity and matter Lagrangians. The dependence on 
the extra dimension and on the details of the gravity and matter models manifests 
itself through the dependence of the function X{z) on the coordinate z, dehned via 
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dz^ = n_(i|/^/(f2+a^(|/)), and on the fnnction RL = a^/a. Focnsing now on Eq. (l5^ and 
redefining X = d ~Y, we find that Eq. fl3^ can be written as 


- n. + v,,f{z)Y = pV, 


(35) 


where 




(36) 


It is now immediate to see that the differential operator on the left-hand side of fl5S]) 
can be written as the prodnct where S = ^ — This decomposition implies 

that S^S is a non-negative operator, ensnring that > 0 and gnaranteeing stability of 
the gravity sector. The zero mode arises for p = 0 and can be readily solved nsing fl33l) 
to obtain 

X{z) = Xo + r .r^dz' , (37) 

J ci[z'y^ 

where Xq and C are integration constants. In order to avoid pathologic behaviors 
as one moves away from the brane, it is necessary to set C —>■ 0, which tnrns X{z) 
into a constant. The constancy of X{z) for p = 0 shows that the zero mode of tensorial 
pertnrbations is the same for all Palatini theories of the form discnssed in this work. This 
pnts forward that the stability of the zero modes observed in the case of GR is a generic 
prediction insensitive to the details of the scalar field and gravity Lagrangian chosen, 
thns snpporting the robnstness of braneworld scenarios against tensorial pertnrbations. 


5. Summary 

In this work we have investigated the backgronnd strnctnre and tensorial pertnrbations 
for generic braneworld scenarios in {d -|- 1)-dimensional geometries with a single extra 
dimension of infinite extent. The family of theories of gravity that we considered is 
constructed trading the Ricci scalar R = g°‘^Rap with a generic function Fi^gap, Rajs) 
and assuming that the metric and affine structures are a priori independent (metric- 
affine or Palatini formalism [9]), an aspect motivated by condensed matter systems with 
defects. Remarkably, this has led to a generic set of ghost-free second-order equations 
for the background, regardless of the particular form of the gravity Lagrangian. We have 
shown that by a suitable choice of variables, the scalar field and background equations 
can be written as a set of first-order equations for the warp factor, H = —W (0)/(d — 1), 
and the scalar field once a superpotential function W (0) is introduced. This property 
suggests that supersymmetric extensions of these theories might be possible. 

We have also shown that the gravity sector in these theories is stable under linear 
perturbations and develops a zero mode (the graviton) bound to the brane with exactly 
the same properties as the massless gravitons of GR. Remarkably, these results are 
robust in the sense that they hold no matter the specific forms of both the gravitational 
Lagrangian, Cq, and scalar field model with generalized kinematics, C{(f),X), chosen. 
The results presented here can be used to study specific models, such as Born-Infeld 
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gravity and its extensions [13 CS]) as well as other theories considered in recent 
investigations |25]. We shall further report on this elsewhere. 
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